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Abstract 

I compute the difference of vector and axial vector current correlators in the weak coupling 
phase of (lattice-regulated) SU(3) gauge theory with two flavors of symmetric-representation dy- 
namical fermions. This is a walking theory at the bare parameter values chosen for the simulation. 
Otherwise, it is not a conventional technicolor candidate. The correlator difference shows scaling 
behavior in the fermion mass, and vanishes in the fermion zero mass limit. Consequences for the 
phenomenology of similar systems which might be candidates for beyond Standard Model physics 
are discussed. I check my methodology against ordinary QCD, by computing the Gasser-Leutwyler 
coefficient Lio and the charged - neutral pion mass difference from an approximate parameterization 
of the correlator. 
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I. INTRODUCTION AND BACKGROUND 



We hope to discover new physics at the Large Hadron Colhder. But even if a direct 
observation of physics beyond the Standard Model is difficult, it still may be possible to 
detect new physics through its effect on Standard Model processes, via precision electroweak 
measurements. One classic observable is the S-parameter of Peskin and Takeuchi[l|. In the 
simple system I will analyze, it is related to the difference between the momentum space 
vector and axial current correlators at zero momentum 



n^.(g) = J d'qexp{tqx){j;:{x)J:^{0)) 



(1) 



The S-parameter is proportional to the limiting value of d{q^'nJ^^{q^)) / dq^ at small g^, after 
Goldstone boson effects are subtracted. 

Recently, many groups have begun to use lattice methods to study candidate beyond- 
Standard Model systems which replace the fundamental Higgs field by new, nonperturbative 
physics [2]. These models may be realizations of technicolor (for reviews, see Ref. jll), or 
they may correspond to "hidden sector" [J] or "unparticle" systems [5] where the new physics 
is approximately conformal. To decide whether any particular model is viable, it is necessary 
to compute observables. This paper describes the measurement of 11^^ for one candidate 
theory, S'f/(3) gauge fields coupled to Nf = 2 fiavors of fermions in the symmetric (sextet) 
representation. Before lattice simulations began, this model was one of many candidate 
theories for walking technicolor j6|-@]. Previous lattice simulations of this system include 



Refs. |10l-ll5j. In this system, Il^,y{q) involves the correlation function of two fermionic (ip) 



bilinears. In my conventions = ip'^^^^—^Tp and = ij^'y ^^"^^j^ip ■ 

The specific lattice calculation is done at a set of bare parameters at which the gauge 
coupling runs slowly or "walks" : the theory is approximately conformal in the zero fermion 
mass limit. Tuning the fermion mass away from zero explicitly breaks conformal symmetry. 
What happens to Uff^lq"^) in this situation is, as far as I know, unexplored. Unfortunately, 
this system is not a candidate for conventional technicolor. At the parameter value where I 
did the simulations, it shows no sign of spontaneous chiral symmetry breaking. 

There is some discussion in the literature of Uff^lq"^) for technicolor candidates which 
have slowly running or "walking" couplings, or which are conformal [isj. In these systems, 
conformal symmetry is is usually broken, either explicitly or through the coupling of the new 
physics sector to the Standard Model. Here, the new physics at very short distance which 
generates a fermion mass, is a bare fermion mass in the lattice action. The question is then, 
what happens to n|^^ as the mass is tuned to zero. The expectation is that n|^^ falls to 
zero in that limit. A second expectation often seen in the Beyond Standard Model literature 
is that can be computed by saturation by a few light resonances in the appropriate 
channels, by tuning their masses and couplings. I will test these expectations. 

By "slowly running" I mean that a suitably defined running coupling constant shows 
small variation, over the range of length scales accessible to a lattice simulation in finite 



volume at a fixed value of its bare couplings. In Ref. [13[, this is done using a lattice version 
of the background field method, in which the system size L represents the scale at which a 
running coupling g'^{L) is measured. (It is called the Schrodinger functional method. Part of 



the extensive literature of the Schrodinger functional include Refs. (17H20||.) It was observed 
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\l3\ that the couphng ran more slowly than perturbation theory predicted. At one loop, this 
is 

1 2b 

— — - = — ^ log sL + constant. (2) 

where bi = 13/3 for this theory. For a scale factor s = 2, this is a change in l/g'^{sL) of 
about 0.038. compared to the measured value at this simulation's bare coupling of l/g'^{L) 
of about 0.38. 

This is not what happens in ordinary QCD. With two flavors of fundamental representa- 
tion fermions, 6i = 29/3 and in perturbation theory the change in the inverse coupling over 
a scale factor of two is about 0.084. In practice, at coupling values used in standard QCD 
simulations, the change is much greater. (Ref. [23| compares their Schrodinger functional 
coupling to perturbation theory.) In fact, the whole framework of high precision lattice QCD 
simulation is built on the assumption that one can do a simulation in which the coupling is 
perturbative at the shortest available distances, and the system will become nonperturba- 
tive at the longest ones. This can can be observed in (for example) a heavy quark potential 
behaving as V{r) ~ 1/r at small r and ar for large r. The ever bigger lattices used in 
state of the art simulations are present both to push the simulation volume (measured in 
centimeters^) ever larger, while simultaneously shrinking the lattice spacing ever smaller. 
"Eliminating lattice artifacts" is a coded phrase for writing down a theory at the cutoff 
scale which is continuum QCD, up to small and calculable corrections, which implies that 
it is weakly interacting there. 

The slow running of the coupling for this model means that on any lattice size accessible 
for numerical simulations, the zero bare quark mass system is conformal for all practical 
purposes. What does control the correlation length (inverse of a mass) is the fermion mass. 
The scaling of the correlation length with fermion mass is given by a critical exponent ym 
which is related to the anomalous dimension 7^ of the operator ipip, defined through the 
running of the fermion mass 

= -lU9')m{f.), (3) 

The quantity 7^ is the interesting parameter for technicolor dynamics. Technicolor sce- 
narios, conventional or not, are said to prefer to have 7™, ~ 1 to simultaneously generate 
phenomenologically interesting fermion masses while suppressing flavor changing neutral 
currents. (Compare the recent discussion in Ref. |2l |.) Our theory is not very desirable 
from a phenomenological point of view: two studies |l2L Il3l | show that 7^ is small, about 
0.35 at the parameter values of the simulation. 

Finally, to complete the list of undesirable features of the system, as one tunes the bare 
gauge coupling larger and larger, it undergoes a first order transition into a confining phase 
where the axial Ward identity quark mass never vanishes. It is not known whether this 
transition is a lattice artifact, or not. For the purpose of this paper, the resolution of this 
question is only of indirect importance. At its bare coupling value, the lattice system has a 
slowly running coupling constant, no chiral symmetry breaking, and no confinement. Let us 
just regard it as a template for some exotic new physics scenario, and see what it produces 
for an electroweak observable: what is n|^^(g^)? 

Readers should note: most lattice calculations are about numbers and precision. This 
calculation has neither. Rather, it is a qualitative study of a theory which does not resemble 
QCD. To try to discover what it does resemble, I (mostly) compute Hj^n rather than the S 
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parameter, because at the bare couplings where I did the simulation, there is no 5* parameter. 
The true electroweak observable S involves subtracting the contribution of the Higgs field - 
or of the particles which replace the Higgs field - from U!^^. There is no evidence of Higgs 
- like dynamics in the weak coupling phase of this theory. 

The realization of Eq. [T]in a lattice simulation involves a long stream of annoying technical 
problems which must be overcome. Fortunately, there are already two lattice calculations 
of the S parameter in QCD, Refs. and [23|, plus earlier lattice work on current-current 



correlators [2J, |25| which provide pretty explicit directions to follow. The problems to be 



addressed are: 

• Lattice artifacts in Eq. [TJ in a lattice calculation 

= P^.m^iQ) + P^.mLiq) + ... (4) 

where Pt{<i) and PlIq) are lattice analogs of the transverse and longitudinal projectors 
and the dots represents additional momentum-dependent terms, proportional to higher 
powers of products of q. 

• If the lattice currents are not conserved, there is an additional quadratic divergence 
S^.uQ/a^ in Il^,u{q)■ 

• If the lattice currents are not local, there are additional contact terms in H^;y(g). 

• There is a lattice - to - continuum regularization factor for each current. If it is different 
for the vector and axial vector currents, the analysis of a quantity depending on their 
difference becomes more fraught. 

If we are only interested in the difference of vector and axial vector currents, the use 
of valence fermions which encode exact chiral symmetry (overlap or domain wall fermions) 
alleviates the second, third, and fourth of these problems. 

The next section describes technical details of the simulations. Readers uninterested in 
them should jump to the following two sections, for qualitative comparisons of Uf^^ from 
ordinary QCD and the walking theory. 



II. NUMERICAL TECHNIQUES AND BACKGROUND 

I will make qualitative comparisons of lattice data from simulations of ordinary QCD 
and of SU{3) gauge theory coupled to two flavors of fermions in the symmetric (sextet) 
representation. In all these studies, the valence Dirac operator is taken to be the overlap 
operator |26l. l27l| . Details of the particular implementation of the action are described 
in Refs. |28l-l32| : suffice it to say that the massless overlap operator is defined as D = 
i?o(l + d{—Ro) I \/ d'^{—Ro)d{—Ro)) where d{m) = d + m for some lattice approximation 
("kernel") d to the massless continuum Dirac operator. The only new ingredient is the 
application to symmetric-representation fermions, already described in Ref . (l2|| . Eigenvalues 
of the squared Hermitian Dirac operator D^D are computed using the "Primme" package of 
McCombs and Stathopoulos[3^^ and are used to precondition the calculation of propagators. 

The lattice analog of Eq. [1] is computed using the difference of improved currents, that 
is, the vector current is 

= - |f-)^2. (5) 
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The axial current is identical, apart from the substitution of 7^75 for 7^. These currents are 
not conserved, and so the correlator of each current is quadratically divergent. However, 
because they are related by a Ward identity, the quadratic divergence cancels in the vector- 
axial difference. In addition, both currents have the same lattice-to-continuum regulator 
renormalization factor (Z-factor). The currents are local enough that U^i, does not have 
contact terms. 

In practice, the correlator is computed using point currents and the "shifted" propagator 

In free field theory and to all orders in perturbation theory 11^^ = at rUg = 0, basically 
because {75,!)"^} = 0. 

This gives a sensible and useful Il^jf - a correlator of local currents with a common Z 
factor, and with the decomposition of Eq. [1] for small momenta. I now have to deal with the 



lattice artifacts in the decomposition of Eq. HI Following Ref . j22| , I observe that they are 
small. This is quantified via the observable 

A.(g) = - J-sHM (7) 



In a slight variation on the method of Ref [22|, I define the appropriate variable as 



(2/a) sing^a/2). (Of course, = {2tx/ L)n^ for integer if the (x direction exhibits periodic 
boundary conditions and its length is L.) A 11^^ which is a superposition of pure longitudinal 
and transverse terms (the continuum decomposition of Eq. [1] with g's replaced by g's) will 
give A = 0. When I come to analysis I will show that A is very small compared to 11^^. 

Next, I have to perform the decomposition of II^i, into YIl and 11^. An easy way to do 
this is to assume that 

^,u{q) = Pl{q)^T{q) + Pi.mL{q) (8) 

where 

P^Aq) = f^t^y - (9) 

and 

Piuiq) = %qu (10) 

are the transverse and longitudinal projectors. Next I form the chi-squared function, in- 
dividual momentum mode by momentum mode, taking the sixteen (/it/) correlators as the 
quantities to be fit, 

x\q) = $^(n^.(g) - Pl{q)TiT{q) - PiAq)^L{q))\ (11) 

and treating 11^ and TIl as fit parameters. Because and are projectors, Tr P^P^ = 0, 
this minimization reduces to simple definitions of weighted averages, 

^ , , j:,.p^um,M 



^L{q) = — 



3(g2)2 

fj,u ( 



(12) 
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Of course, all lattice data from the same set of configurations are highly correlated. The 
uncertainties in the pictures which follow are computed by folding Eq. [12] into a single 
elimination jackknife: I delete a lattice from my data set, compute 11^, Ul, and A, and then 
present a jackknife average and uncertainty. 



III. ILLUSTRATIONS FROM QCD 

The remainder of the paper is a qualitative comparison of 11^^ from ordinary QCD and 
from sextet QCD. I begin with ordinary QCD. 

I have two data sets. The first is a quenched set of 20 lattices computed using the overlap 
operator for valence quarks on a background of Wilson gauge action configurations at a 
gauge coupling (3 = 5.9. The lattice size is 16"^ sites. The lattice spacing is about 0.13 fm 
from the rho mass or 0.11 fm from the Sommer parameter. This data set can be compared 
to an extensively- analyzed companion used in a calculation of Bk, the kaon B-parameter 
'3^ . (The parameters of the kernel action used in this study differ by about one per cent 



from those of the Bk project, but that is not going to matter given the qualitative nature 
of my presentation.) 

The second data set is a smaller volume (12^ sites) set of simulations with Nf = 2 flavors 
of dynamical overlap fermions. The data set had three quark masses, arriq = 0.03, 0.05, 
0.10, at a lattice spacing of roughly 0.14 fm from the Sommer parameter. These lattices 
have been used in several small projects [35!] to date. They are not really big enough for 
reliable spectroscopy measurements. 

I cannot extract an 5* parameter from these results (though I tried). The reason is that 
to get S from Ilf^^ requires doing a fit of the data to a chiral perturbation theory calculation 
and removing the single pseudoscalar propagator contribution. S is then related to one 
of the Gasser-Leutwyler O(p^) coefficients. Doing this requires going to very low g^, so 
that chiral perturbation theory is in its domain of applicability. This in turn requires with 
present techniques making the lattice volume large. (For example, Ref. used a 16'^ x 32 



lattice.) My 16^* lattice is large but it is a quenched data set and I am unaware of the 
appropriate quenched chiral perturbation theory calculation. And of course, the quenched 
approximation is obsolete as a potential high-precision venue. For this study, I do not need 
S. 

Figs. [1] and |2] show the A parameter and Uff^ from the quenched data set. In keeping 



with past practice [22|, |23|, I have plotted g^II;^^. A seems satisfactorily small compared to 
n^^. Unsurprisingly, 11^^ qualitatively resembles other published results. 

The superimposed curves in Fig. |2] show phenomenological parameterizations of the data, 
where II^''^ is saturated with a set of resonances. Introducing two quark flavor labels to 
characterize the currents, the pseudoscalar decay constant is defined through 

(0|M7o75C?|vr) = m^f„ (13) 
(so /tt ~ 132 MeV) while the vector meson decay constant of state V is defined as 

{0\u^id\V)=mlfve, (14) 
and the axial vector meson decay constant of state A is 

{0\uja5d\A) = mlfAei. (15) 
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FIG. 1: (f'/Sy~^ vs (f' vs (f' from quenched overlap fermions at /3 = 5.9. Valence masses are 
(octagons), aniq = 0.10 (squares), arriq = 0.05 (diamonds), arriq = 0.035. (crosses) arriq = 0.025, 
and (fancy diamonds) aruq = 0.02. 




FIG. 2: q^I]}^~^ vs from quenched overlap fermions at /3 = 5.9, As in Fig. [U valence masses are 
shown by octagons for amq = 0.10, squares for arUq = 0.05, diamonds for arUq = 0.035, crosses for 
arUq = 0.025, and fancy diamonds for arUq = 0.02. The lines are the model parameterizations. 
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FIG. 3: As in Fig. [2l but a logarithmic scale to show the approach to Weinberg sum rules. 



ej is a unit polarization vector. The transverse current correlator is then 



In a theory with spontaneous chiral symmetry breaking, the couplings are constrained by 
the first 

Y.fv^^v'Y.fAM\~fl = Q (17) 

V A 

and second 

E/^<-E/l^A = o (18) 



Weinberg sum rules 36|. When the two sum rules are satisfied, I\.J^^{q^) falls off asymptot- 
ically in q as 

A common phenomenological approach is to saturate the correlator by the lowest states 
in the vector and axial channels, the pion, rho and ai meson[3T]. This involves fixing five 
parameters, the three couplings /^r, /p, and fa^ and two masses nip and nia^. In principle, 
all of these can be measured in lattice simulations. In practice, fa^ and involve difficult 
measurements because the signal to noise ratio in a mesonic correlator (whose mass is M) 
scales like exp((M — 2mT,)t). 

I attempted to extract all five quantities from spectroscopic fits to the quenched data. 
I could do this only for the heavier masses in the data set. Results are shown in Table 
[H Using them in the correlator does not reproduce either the data or the Weinberg sum 
rules. (Perhaps this result is neither surprising nor controversial.) However, I do want to 
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mq 


fn mp 


fp 


mai fai 


fit fit /ai 


Iio X 10^ 


a'^Aml X 10^ 


0.020 


0.085(1) 0.55(2) 


0.282(5) 




0.657(15) 0.197(3) 


5.1(6) 


0.62(4) 


0.025 


0.089(1) 0.56(2) 


0.277(4) 




0.684(14) 0.185(9) 


5.3(6) 


0.66(5) 


0.035 


0.095(1) 0.57(1) 


0.271(3) 




0.723(8) 0.168(6) 


5.6(7) 


0.71(3) 


0.050 


0.100(1) 0.60(1) 


0.261(2) 


0.93(5) 0.182(3) 


0.780(8) 0.154(5) 


5.5(8) 


0.81(3) 


0.100 


0.1168(5) 0.689(4) 0.235(2) 


1.00(3) 0.170(2) 


0.995(10) 0.113(4) 


5.3(6) 


1.16(2) 



TABLE I: Measured lattice parameters (in lattice units; everything is scaled by a factor of the lattice 
spacing a and the overall Z factor is left out) from quenched 16^ simulations, used in constructing 
the parameterization of the transverse correlator. The last four columns are fit results from forcing 
the spectral function to obey the first and second Weinberg sum rules, saturated by the rho and 
ai mesons. 



ruq 


/tt mp fp 


fit fit /ai 


Lio X 10^ 


a'^Aml X 10^ 


0.03 
0.05 
0.10 


0.092(3) 0.57(2) 0.22(2) 
0.135(5) 0.64(2) 0.25(2) 
0.138(6) 0.80(1) 0.23(2) 


0.84(10) 0.10(3) 
1.19(26) 0.07(4) 
1.20(15) 0.10(3) 


4.8(4) 
7.1(7) 
5.3(5) 


0.81(12) 
1.24(23) 
1.64(18) 



TABLE IL Measured lattice parameters from dynamical 12^ simulations, used in constructing the 
parameterization of the transverse correlator. The last four columns are fit results from forcing 
the spectral function to obey the first and second Weinberg sum rules, saturated by the rho and 
ai mesons. 



compare the QCD results to those from the walking theory. Accordingly, one approach to 
parameterizing the data consists of taking the values of rrip, and fp from fits to lattice 
data, and determining /^^ and rria^ by forcing a solution to the to Weinberg sum rules. 

The result of this calculation is shown as the curves in Fig. [2l In this figure, the bursts at 
q = are just the values of Fig. |3]shows g^II^"^ on a log-log scale to expose the approach 
to scaling at large q. The data qualitatively resembles the analytic parameterization. The 
Weinberg sum rules are statements about the asymptotic behavior of the current correlators. 
Lattice QCD correlators presumably attempt to encode these statements until they are 
overwhelmed by lattice artifacts. 

Next, we turn to figures from the 12'* volume dynamical overlap: Fig. H] shows the A 
observable. Fig. [5] shows q'^Il^~^ vs q^. Again, the superimposed curves are from Eq. [161 
taking fit values of f-,,, fp and nip, and determining fa^ and m^^ by saturating the two 
Weinberg sum rules. Fig. |6] reproduces Fig. [5] on a log-log scale to show the Weinberg sum 
rules at work. As a qualitative parameterization of the data, this procedure works well. 

I cannot resist using the pseudoscalar and vector decay constants, plus the vector meson 
mass, to compute the Gasser-Leutwyler parameter Lio and the charged - neutral pion mass 
difference Am^. via the Das, Guralnik, Mathur, Low, Young sum rule jssf. In the "lowest 
mass dominance" approximation to the spectral function, Lio = {fp — = fp{R^ — 1) 

and 

A 2 3a 1 , m2 3a 

^"^^ = 4^ (^-^) ^ = ^ 

where R = 1 ~ {fn/i'm.pfp))'^ is given to remind the reader, what are the real independent 
variables. The results are shown in the two tables and in Fig. [71 A fit to a linear dependence 
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FIG. 4: q^A^ ^ vs vs from dynamical Nf = 2 simulations. Valence masses are (octagons) 
aniq = 0.10, (squares) aniq = 0.05, (diamonds) aruq = 0.03. 



0.005 



0.000 



-0.005 



-0.010 



-0.015 



0.020 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




















1 1 1 1 


1 1 1 1 1 







FIG. 5: g^n^ ^ vs from dynamical Nf = 2 simulations. Valence masses are (octagons) arUq 
0.10, (squares) arUq = 0.05, (diamonds) arUq = 0.03. The lines are the model parameterizations. 
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FIG. 6: As in Fig. [5l but a log scale to show the approach to Weinberg sum rules. 



on the quark mass gives Lio = 5.1(6) x 10~^ and 5.3(5) x 10~^ for the dynamical a nd q uenched 
data sets. This is in good agreement with large scale simulation results of Refs. U 23 1 and 
with phenomenological estimates (39| . 

Fitting the dynamical data set's Am^ in the same way gives aAm^ = 0.49(18) x 10^^ 
or (with a = 0.14 fm) 980(360) MeV^. The quenched data has aAm^ = 0.48(3) x lO^^ for 
1540(100) MeV^ or 1100(70) MeV^ taking a = 0.11 or 0.13 fm. These results agree nicely 
with experiment, 1261 MeV^. Readers should note that while these are lattice results, they 
are indirect and depend on the additional assumption that the spectral functions can be 
saturated by three input parameters plus the two Weinberg sum rules. 

The summary of this section is that both quenched and dynamical lattice QCD possesses 
a g^n^~"^(g^) which at small q remains finite as the quark mass vanishes (it reduces to 
— /^). At large q it vanishes roughly in accordance with the Weinberg sum rules. A simple 
three-parameter function using information from the lowest excitations in the axial and 
vector channels seems to reproduce the data reasonably well and produce phenomenologically 
sensible results. 



IV. FROM SEXTET FERMIONS 



I repeat the calculation of g^n^~'^(g^) using sextet representation overlap fermions on 
a background of dynamical simulations of Nf = 2 flavors of clover fermions in the sextet 
representation and SU{3) gauge fields. The dynamical simulations were performed using 
the Wilson gauge action at a coupling /3 = 5.2 and a hopping parameter k = 0.1285. This 
corresponds to an AWl (axial Ward identity) quark mass of am^ = 0.044. 1 used five 
valence masses, am„ = 0.100, 0.050, 0.035, 0.025 and 0.020. The lattice volume is 16"^ sites. 
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FIG. 7: Panel (a) shows lattice data for Lio, and panel (b), for Am^ from the analysis described 
in the text. Squares show quenched data and crosses are for Nf = 2 simulations. 

At this set of simulation parameters, the system seems to be deconfined (the string tension is 
immeasurably small) and chiral symmetry is restored (observed through the parity doubling 
of states). 

Fig. [8] shows the A observable. It seems to be acceptably small. Fig. [9] shows g^H^""^ vs 
g^. Fig. [TOl shows —(fY\^~^ vs g^, replotted with logarithmic axes. These figures show very 
different behavior from the case of QCD. 

Notice that in the large-g^ limit, 11^^ still falls to zero. This is consistent with the 
Weinberg sum rules, and is expected for any theory with an ultraviolet-attractive fixed 



point at q = |40 



Just as in the case of QCD, it is possible to measure spectroscopy in the pseudoscalar, 
vector, and axial vector channels. It is also possible to measure the three decay constants. 
In fact, it is easier to extract the axial vector meson properties here, than it is in an ordi- 
nary QCD simulation. This is because the system is chirally restored and the pion is not 
appreciably lighter than the other states. (The same effect occurs in the high temperature 
phase of QCD.) These parameters are shown in Table lllli The obvious feature they show 
is chiral symmetry restoration: note the near degeneracy of vector and axial vector masses 
(called TUp and in hadronic analogy), which becomes more pronounced as the fermion 
mass vanishes) and of the vector and axial vector couplings. 

A quick check with a calculator shows that the three lowest states with their couplings 
do not saturate either Weinberg sum rule. So 1 play the same game that 1 did for the QCD 
data sets, determining rria^ and fa^ by forcing the two sum rules to saturate. The result of 
this exercise is shown in Fig. [10] (for the = 0.100 data set, only). This fails completely. 

1 was somewhat surprised by this result, since lowest mass dominance is a venerable 
ingredient of strongly coupled beyond - Standard Model phenomenology (compare Ref. [H, 
or for a recent review, Ref. (4l|), even sometimes of conformal beyond - Standard Model 



phenomenology (compare Refs. [42, |43|. However, note that 1 am looking at a system in 



which the breaking of (near) conformal symmetry is by a small mass, and 1 am looking at 
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FIG. 8: q'^A^~^ vs from valence sextet overlap fermions on background sextet dynamical 
clover simulations at /3 = 5.2, k, = 0.1285. Valence masses are (octagons) aniq = 0.10, (bursts) 
aniq = 0.075, (squares) aruq = 0.05, (diamonds) aniq = 0.035. (crosses) aruq = 0.020. 
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FIG. 9: q^Il^~^ vs from valence sextet overlap fermions on background sextet dynamical 
clover simulations at /3 = 5.2, k. = 0.1285. Valence masses are (octagons) aniq = 0.10, (bursts) 
aniq = 0.075, (squares) aniq = 0.05, (diamonds) aniq = 0.035. (crosses) aniq = 0.020. 
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FIG. 10: As in Fig. [U but a log- log scale. The line shows an attempted parameterization of the 
data, forcing the two Weinberg sum rules to be saturated by the three lowest resonances. 



0.020 
0.035 
0.050 
0.075 
0.100 



m, 



■ax 



fai 



0.00050(4) 0.33(1) 0.0046(2) 0.38(1) 0.0044(2) 
0.00080(6) 0.37(1) 0.0045(2) 0.42(2) 0.0043(1) 
0.00102(8) 0.37(1) 0.0044(2) 0. 47(2) 0.0042(2) 
0.00131(9) 0.43(1) 0.0041(1) 0.54(2) 0.0037(10) 
0.00151(10) 0.48(1) 0.0037(1) 0.61(2) 0.0037(2) 



TABLE III: Measured lattice parameters from valence overlap fermions from sextet QCD simula- 
tions. 



rn^/q^ small. Often, the phenomenology is concerned with the other limit, im?/q^ large. 
I convert n^""^ into an "effective" S-parameter by numerically differentiating it. 



1677 



(20) 



(For the resonance model, Eq.[Tni this would give S^ff = 167r(/y — /^).) This I do in the most 
naive way, sorting the data into increasing and just taking successive differences. This 
rapidly becomes noisy at bigger g^, so I arbitrarily kept the ten smallest values, q < 1.24, 
for each rriq. I plot the result in Fig. [TTJ The collapse of the data to zero with the fermion 
mass is obvious. Presumably, a true "new physics" prediction would multiply Seff by an 
overall scale representing the coupling constant of the dynamics to the W or Z boson. 

It is a worthwhile exercise to look for scaling behavior in the data. To do this, I replot 
the data in Fig. [T2i The data seems to fall on a scaling curve in terms of the dimensionless 
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FIG. 11: Derivative of the sextet QCD correlator Seff = 167rA(g^n^~"^)/Aq^.. Valence masses 
are (octagons) arug = 0.10, (bursts) aniq = 0.075, (squares) aniq = 0.05, (diamonds) anig = 0.035. 
(crosses) anig = 0.020. 



ratio q /m . In free field theory, S^ff = (48/7r )(m /g ), which seems to be a reasonable 



parameterization of the data. However, motivated by various didactic reviews of the 
expected algebraic scaling behavior of conformal theories, I also tried a little fit to 

2 

77? 

Seff = a{^r. (21) 

The line is the result, p = 0.953(3). The uncorrelated chi-squared is 924 for 38 degrees of 
freedom, but I should point out that the data is in fact strongly correlated. 

Again, I remark that this is not the usual case studied in the technicolor literature, where 
the fermion mass is often taken to be large compared to q. Presumably to apply this result 
to some hidden sector beyond - Standard Model phenomenology, one might take q^ = m|. 
(Or for another alternative, see Ref. (isl].) 

I believe that there are lessons to be drawn from this exercise. First, because the system 
only has explicit chiral symmetry breaking at nonzero fermion mass, Il^~^{q'^) does vanish 
at zero fermion mass. Il^~^{q'^) falls to zero for any value of the mass, at any nonzero q. 

Second, a simple parameterization of Il^~^{q'^) in terms of the (measured) properties of 
the lowest resonances fails, even though away from rUq = the theory is not conformal; 
it is a theory of resonances. Presumably, one must just include more of the excited state 
spectrum in the sum. S^ff shows power law scaling in terms of rn? /q^. 
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FIG. 12: Derivative of the sextet QCD correlator Seff = l6TTA{q^U}^~^)/Aq'^ plotted as a function 
of the scaling combination if' jrr? . Valence masses are (octagons) am^ = 0.10, (bursts) arriq = 
0.075, (squares) aruq = 0.05, (diamonds) arug = 0.035. (crosses) arug = 0.020. The hue is the 
straight power law fit, Eq. [2TJ 



V. CONCLUSIONS 

This paper was a qualitative investigation of a precision electroweak observable for a 
lattice model of gauge fields and two flavors of symmetric representation fermions. Before 
lattice simulations began, this model was one of many candidate theories for walking tech- 
nicolor 0-[9|. At the bare parameters where 1 did my simulations, it does not present the 
appearance of being a conventional walking technicolor model. However, it does show con- 
formal behavior "for all practical purposes," that is, at the parameters of the simulation the 
scale dependent gauge coupling runs very slowly [l3|. 

1 only asked simple questions, because the literature for these models is relatively sparse. 
However, 1 discovered (as expected) that H^~"^(g^) vanished for all q as the quark fermion 
mass vanished, and that it vanished at large q for all values of the fermion mass. The 
effective S-parameter is a scaling function of m?/q^. 1 was surprised to discover that a 
phenomenologically popular parameterization of H^~"^(g^) in terms of the properties of the 
lowest resonances did not reproduce the data. 

It might be, that all models are just different, and what I found for H^~"^(g^) is not 
relevant to any other system. However, there are other would-be technicolor candidates in 
the literature, which share some features with this model for regions of their bare parameter 
spaces. n|^^ in "minimal walking technicolor" (gauge group SU (2) with a doublet of adjoint 
representation fermions j46|) is an interesting target for a lattice study. 

And of course, phenomenologists should expect more calculations of H^~^(g^), done by 
groups on models which are (hopefully) more attractive technicolor candidates - perhaps 
ones which are not so conformal as this one. 
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